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ABSTRACT 

Using our extensive R OSAT X-ray observations of M33, we confirm a 3.45 day eclipse 
period for the Einstein source X-7 ( Larson fc Schulman 1997 ) and discover evidence for 
a 0.31-s pulse period. The orbital period, pulse period and observed X-ray luminosity 
are remarkably similar to SMC X-l. We therefore suggest M33 X-7 is a neutron star 
high mass X-ray binary with a 15-40 M O/B companion and a binary separation of 
25-33R.0 if the companion is almost filling its Roche lobe. 
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1 INTRODUCTION 

M33 is an Sd spiral galaxy in th e Local Group, which at 
~ 800kpc (van den Bergh 1991), is slightly more distant 
than M31. At an inclination of i ~ 56 °, the disk of M33 
is accessible to observation without a large amount of ab- 
sorption. The first X-ray studies of M33 made with the Ein- 
stein Observatory Imaging Proportio nal Counter ( [ Long et 
al. 198l|) and High Resolution Imag er ( |Market fc Rallis 1983[ 

Trinchieri, Fabbiano & Peres 1988), resulted in the detection 

l 



of 13 sources brighter than 10 ergs - , one of which, X-7, 
was flagged as variable. As a result of their luminosities the 
X-ray sources in M33 are thought mostly to to be super- 
nova remnants or X-ray binaries , and X-7 clearly falls into 
the latter category. Peres et al. (1989) and later, Schulman 



et al. (1993; 1994) discovered eclipse episodes in X-7 and 



proposed a period of 1.78 days. This strongly suggests that 
X-7 is a High Mass X-ray Binary (HMXB). A recent reanal- 
ysis of archival Einstein, R OS A T and ASCA observations of 
X-7 by Larson and Schulman ( 1997 ) argues in favour of a 
3.45 day period. 

Here we discuss a considerably larger set of ROSAT 
observations of X-7 which extend over a time base of 6 years. 
We use these observations to confirm the 3.45 day period and 
report the discovery of a 0.31 s pulse period, which allow us 
to establish the X-ray source as a pulsating neutron star. We 
have used this dataset previously to d iscuss the tempora l 
variability of the nuclear source X-8 (Dubus et al. 1997) 



(hereafter paper I). A more general discussion of all of the 
sources will be the topic of a future publication. 



2 OBSERVATIONS 



Multiple observations using either the ROSAT High Res- 
olution Imager (HRI) or Position Sensitive Proportional 
Counter (PSP C) w ere carried out between 1991 and 1997 
(see Triimper (L984) for instrumental details). Most of these 
were centered on the optical nucleus of M33 at a l h 33 m 50!4, 
£ 30° 39' 36" (J20 00). With a luminosity of » 10 38 erg/s 
([Long et al. 1996; ), X-7 is easily detected in all cases. The 



different observations are presented in table [l]. The total ex- 
posure time is 233-ks of HRI (PSPC: 50-ks) nucleus pointed 
observations and about 92-ks of off-centered HRI data. From 
the long set of HRI observations pointed at the nucleus of 
M33, we get a best position for X-7 of q(J2000) l h 33 m 34!2 
and 5 30°32'08'/7, with an error of 5" to 10". This is based 
on deriving a correction for the ROSAT X-ray position of 
t he nucleus by using the accur ately known optical position 
( |de Vaucouleurs fe Leach 1981 ). 



The barycentre-corrected photon arrival times were ex- 
tracted from a 1' region centered on the best position for 
X-7. Background counts were also extracted from a 2' ra- 
dius disk centered on a l h 33 m 20!7 and 5 30°31'5l'/3 (J2000), 
close to X-7 and at roughly the same angle from the axis. 
This choice was dictated by the presence of another source 
close enough to X-7 that a standard background annu- 
lus would have been unreliable. Depending on the analy- 
sis method, we used the data both in arrival time format 
and in spacecraft orbital averages. For the latter, we group 
together good time intervals within 3-ks of each other (~ 
orbital time within which M33 is in the field of view) and 
correct for b ackground and vignetting as described in David 
et al. ( 1996). After corrections, the mean HRI count rates 
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Table 1. ROSAT observations of M33 X-7 



Obs. 


Offset (') 


Dates 


Duration (ks) 


rp23a 


8 


1991 29-30 Jul 


29.1 


rp23b 


8 


1992 10 Aug 


5.0 


rp2Sc 


8 


1993 7-9 Jan 


16.3 


rh20a 


7 


1992 8-12 Jan 


19.1 


rh20b 


7 


1992 1-3 Aug 


15.8 


rh86 


12 


1994 1-11 Aug 


16.4 


rh87 


8 


1994 27 Jul - 7 Aug 


30.5 


rh88 


6 


1994 27 Jul - 7 Aug 


24.5 


rh89 


8 


1994 6-8 Aug 


20.3 


rhGOa 


8 


1994 10-11 Aug 


8.0 


rh46 


8 


1995 15-15 Jan 


24.6 


rh60b 


8 


1995 10-16 Jul 


40.9 


rhlln 


8 


1996 18 Jan - 8 Feb 


46.4 


rhlla 


8 


1996 17-27 Jul 


44.6 


rh03 


8 


1997 10-14 Jan 


33.9 



The first three observations were made by the PSPC while the 
others were carried out with the HRI. Offset is the angle between 
X-7 and the pointing center (the nucleus of M33 in most cases). 
The abbreviations used correspond to the last digits of the actual 
sequence number. 

were about 11 counts ks -1 (PSPC, 39 counts ks - ) for X- 
7 and about 3.5 counts ks -1 (PSPC, 1 count ks -1 ) for the 
background. Deadtime corrections are negligible at such low 
count rates. 
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Figure 1. Folded lightcurve of X-7 over the best period of 3.45d 
(two cycles are represented). Each point represents an orbital flux 
average. The continuous line is the best fit model. Phase corre- 
sponds to the mid-eclipse time derived in the text. 

we derive the following ephemeris for the mid-eclipse time: 
JD2448631.5 ± 0.1 + 7V(3.4535 ± 0.0005) 
The folded data with the model superimposed is shown 
in figure The mean X-7 count rate in the constant part of 
the model is found to be « 13 counts ks - . 



3 ANALYSIS AND RESULTS 

The analysis methods are similar to the ones applied to the 
nuclear source M33 X-8 described in paper I. We searched for 
power at frequencies between 10 -3 and 10 Hz using a modi- 
fied Fourier spectrum (see Appendix A). On corrected mean 
orbital fluxes, we looked for periodicities using a Lomb- 
Scargle normalised periodogram. 



3.1 The 3.45 day eclipse period 

All of the observations of X-7 which included an eclipse fea- 
ture were flagged as variable with better than 99.9% confi- 
dence. We then removed the data taken during the eclipse 
times and found no evidence for variability, i.e., out of eclipse 
the source appears to be steady. We modeled the folded 
lightcurve from X-7 as a constant flux + linear ingress and 
egress + an eclipse interval with zero flux. Minimization of 
the five parameters (period P, mid- phase of the eclipse m id, 
eclipse duration A^ cclipsc , ingress duration A0 ingrcss , egress 
duration A0 eL 
sponds to JD 2448628.5) 

P = 298.38 ± 0.02 ks 
0mid = 0.88 ± 0.02 A<j!> cclipsc = 0.20 ± 0.03 
A0 ingress = 0.10 ± 0.05 A0 cgrcss = 0.01 ± 0.01 
The errors are estimated from all the minima found 
with xS~l using different minimization routines, initial val- 
ues and number of free parameters. The period of 3.45 
days is in agreement with that found by Larson et al. 
( 1997 ). The values for P, <^ mid and A^ eclipse are quite ro- 



rcss) l e< i to tne following best fit (<f>—0 corre- 



bust and accurately fitted. eas was always close to 

(symetric profiles led to a reduced xt~2) while A0 ingress 
was always >0. Using P, ci m id and the associated errors, 



3.2 Evidence for a 0.31s pulsation 

To search for periodicities at high frequencies, we split the 
PSPC and HRI data into 500s intervals of continuous data 
where X-7 is positively detected (i.e. we exclude the eclipse 
intervals) and calculated the summed Rayleigh power spec- 
trum as described in the appendix. By doing so we lose about 
half of the available data to poor statistics. Smaller inter- 
vals would have resulted in less data being used because of 
the necessity to have at least 5 counts in each interval (see 
appendix A) while there are few continuous data intervals 
longer than 1000s. Choosing an interval length of 1000s in- 
stead of 500s did not substantially change the conclusions. 

The resulting summed power spectrum is binned (re- 
ducing our frequency resolution) to increase signal-to-noise. 
Assuming a \ 2 distribution for the power (see appendix A), 
logarithmic bins in frequency showed no deviations from 
Poisson statistics such as red noise in our range of 10 -3 - 
10 Hz (photons are time-tagged to a resolution of ~10 _3 s). 
However, the data showed a significant (99.9% confidence) 
power enhancement at 3.17±0.08 Hz with linear binning, 
the error being the half-width of the peak in figure 2. The 
signal is washed out with logarithmic binning because of the 
corresponding poorer frequency resolution. Bin sizes of 70 to 
120 independent Fourier spacings gave similar results. 

As discussed in the appendix, we relied on 1000 sim- 
ulated Poisson fluxes to fully assess the significance of this 
peak. The total number of counts n in each simulated flux is 
drawn from a Poisson distribution with the X-7 mean flux. 
Using the properties of Poisson fluxes, the n arrival times are 
uniformly drawn between and the exposure time, sorted 
and then shifted to account for the observation time window. 
The simulated data is processed as with X-7. A Kolmogorov- 
Smirnov test comparing the power distribution of the X-7 
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4 DISCUSSION 

4.1 The orbital lightcurve 

Our well-sampled orbital lightcurve (figure [l]) shows clearly 
that the X-7 eclipse ingress is gradual, but the egress is fast. 
This is compatible with the standard picture of an extended 
region of the disc which precedes the secondary in phase, 
analogous with features well-known in low-mass X-ray bina- 
ries (see e.g. White et al 1995). In high inclination systems 
this extended disc region (created by the impact of the mass 
transfer stream on the disc) will start to obscure the com- 
pact object before its eclipse by the secondary. We therefore 
believe that Roche-lobe overflow plays a significant role in 
X-7, even if it is a HMXB. 



Figure 2. Fourier power spectrum for X-7. The power spectra of 
500-s data intervals from the PSPC and HPJ have been summed 
together. The total power spectrum was then smoothed by av- 
eraging the nearest 100 contiguous frequencies (there are 5000 
independent frequencies in the non-averaged power spectrum). 
The dotted line is the theoretical 99.9% confidence level assum- 
ing a x 2 statistic for the power distribution (the dashed line at 
y=2 being the expected average) and the dots represent the max- 
imum power found at each frequency in the 1000 simulated fluxes 
(i.e. the calculated 99.9% confidence level see section 3.2). 



spectrum with that calculated from the simulations showed 
no deviations, meaning the simulations adequately repre- 
sented the flux from X-7. The maximum power at each fre- 
quency from the 1000 simulations gives an estimate of the 
99.9% confidence level. As expected (see appendix), the peak 
in figure 2 is clearly above the maximum powers found in 
the simulations. This result is unchanged if Poisson fluxes 
are simulated with exactly the same number of counts as X-7 
in each continuous interval (to account for the eclipses). 

Though the feature seems broader than what would be 
expected from a simple sinusoidal pulse, it is not broad 
enough to qualify as a quasi-periodic oscillation (see e.g. 
van der Klis, 1989) and we therefore interpret it as an X-ray 
pulse. Knowing from the lightcurve the eclipse phase and or- 
bital period, we tried correcting the photon times for their 
travel time across the system assuming a circular orbit. This 
did not lead to conclusive results. The Doppler effect and a 
non-zero pulse derivative are not enough to explain the ob- 
served peak width ~0.1 Hz. Assuming a ~50 light-second 
binary separation by analogy with SMC X-l (see discus- 
sion), the Doppler shift is only of the order of 10 -3 Hz. A 
frequency drift of the order of 0.1 Hz over ~5 years implies 
v ~ 10~ 10 Hz s _1 , at least an order of magnitude higher 
than what is observed in galactic X-ray pulsars. 

Our data spans ~5 years, implying a folding frequency 
resolution of ~ 10~ 9 Hz. As the count rate is not high enough 
to detect the pulsation within each observation we cannot 
fold the data and give a pulse shape with confidence. This is 
not in contradiction with the power spectrum analysis which 
does not include phase information (i.e. the distribution of 
observations in time). At this point it is not possible to check 
if the variability is coherent. Although pulsed emission seems 
a reasonable assumption, the power excess could arise from 
variability of a different nature (e.g. broadband variability 
increasing the chance of spurious detection). 



4.2 System parameters 

As a 3.45-d eclipsing system with a 0.31-s pulse period, X- 
7 is very similar to the high mass X-ray binary SMC X-l 
which has an orbital period of 3.89-d, a pulse p eriod of 0.71 
s and an X-ray lum inosity of 5 • 10' 
fc McClintock 1995|). The 10 38 



38 erg s" 1 (van Paradijs 
erg s _1 luminosity of X-7 is 
compatible with the Eddington luminosity of a neutron star 
accreting at ~lO~ 8 M yr~\ As is the case with SMC X-l, 
the companion is likely to be an O or B star, as suggested by 
the loc ation of X-7's error box in the dense OB association 
HS13 (Humphreys & Sandage 198C ). Roche lobe overflow 



is also significant in SMC X-l (van Paradijs & McClintock 



1995), and so we shall assume in our calculations that the 



mass-losing star in X-7 either fills or is close to filling its 
Roche lobe. 

Under this assumption, we take -RrocIic = R.2 to be the 
radius of the star, and adopt a binary inclination of 90° . We 
can now use our observed eclipse duration of 0.20 in phase 
to give R.2 /a — 0.59, where a is the binary separation. Now 
the ratio of the Roche lobe radius to a is (for circular orbits) 
a function only of the mass ratio q [=Mx j 'M%) of the binary 
(Egglcton 1983), from which we derive q = 0.085. Combin- 
ing this with Kepler's third law and an assumed compact 
object mass of 1.4M© (a canonical neutron star mass for 
our observed X-ray pulsar) gives a = 25Rq, M2 = 16.4M© 
and R2 = 14.8Rq. This corresponds to an early-type star, 
and is significantly larger than a ZAMS star (which would 
have a radius of only 5.8Rq for this mass). These parameters 
are consistent with a B0III star which has My = -5.1. 

However, the inclination could be as low as 70°, for 
which our data would then give q = 0.036, a — 32.9R©, 
M2 = 39Mq and R2 = 21R0. Such a star is closer to a some- 
what brighter (My = -6.5) 071 spectral type, and gives an 
indication of the uncertainties involved in this calculation. 
Interestingly, the brightest star in HS13 has V _ 18.8 and, 
with a distance modulus of 24.65 (Humphreys & Sandage 



198C), this corresponds to an absolute magnitude mid- way 



between the B0III and 071 spectral types. 



5 CONCLUSION 

We confirm M33 X-7 has an eclipse period of 3.45 days as 
Larson & Schulman (1997) had found using data from dif- 
ferent satellites. Using a Fourier power spectrum method 
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adapted to low count rates, we also find evidence for a 0.31- 
s period, which identifies the compact object as a neutron 
star. From our best fit to the orbital lightcurve, we find the 
companion is an early type O or B star with M > 15M0 and 
which is likely to be losing mass predominantly by Roche 
lobe overflow. The optical location of X-7 within a dense O- 
B association is consistent with the expected massive com- 
panion. When identified, the optical counter part is likely to 
show e llipsoidal a nd /or heating variations ( |van Paradijs & 
McChn fgck 1995| ). 
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APPENDIX A: TIME ANALYSIS AT HIGH 
FREQUENCIES 

Al Motivation 

Searching for periodicities in our observations with the 
Lomb-Scargle periodogram proved impractical at high fre- 
quencies > 0.001 Hz (the timescale of an orbit). First, the 
periodogram only extends to the Nyquist frequency of the 
data which is inversely proportional to the bin size chosen for 
the time series. The flux from X-7 being rather weak, small 
bin sizes imply long time series of zeros with a few l's at the 
binned photon arrival times. This can be computationally 
problematic to handle. Second, the folding frequency reso- 
lution of ~ 10" 9 Hz for our dataset leads to an enormous 
number of frequencies to search above 0.001 Hz. 

In a Fourier analysis, it is thus customary (e.g. van der 
Klis, 1989) to split long observations into S smaller continu- 
ous intervals of equal length, bin the photon arrival times on 
some multiple of t he satellite clock and then sum the result- 
ing renormalised ( Leahy et al. 1983) Fast Fourier Transforms 
(FFT). With the Leahy renormalisation, the power density 
probability follows a xls l aw when the number of counts 
in each interval is sufficiently large. Splitting observations 
increases the signal-to-noise ratio as compared to a single 
interval power spectrum by a factor vS, though at the price 
of frequency resolution. 
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A2 Rayleigh power 

Binning makes no sense with a small number of events so, in- 
stead of the Leahy power spectrum, we calculated a Rayleigh 
power spectrum for each interval : 

^F n {v)= ^^cos(27ri/t fc )^ + ^^sin(27r^ fc )^ 

if n photons were recorded at times tk in the interval [0, T] 
and v is one of the independent frequencies {1/T, 2/T, ...}. 
The Rayleigh power is essentially the Fourier power of the 
photon arrival times modeled as diracs at tk (infinitely small 
bin size). The only limits to the highest sampled frequency 
are the time-tagging resolution of the instrument or the com- 
putation time requirements. As above, the power density 
probability follows a x\ l aw when the number of counts in 
the interval is sufficiently large (n <; 100; see Brazier, 1994, 
and references therein). 

However, the probability law is already quite close to a 
X% when n <; 5 though the tail of the distributions statis- 
tically differ. In practice, simulations show that the proba- 
bility to have a large value for F n (u) when 5 ^ n S~ 100 is 
smaller than the corresponding probability for a x\ distri- 
bution i.e. the significance of a high peak value in the power 
spectrum will be under-estimated if it is calculated with a 
x\ distribution. If one assumes the x\ i s still an adequate 
representation of the power probability law even with low 
n, then the techniques described by van der Klis (1989) for 
the Leahy renormalised power spectra can be similarly ap- 
plied (summed intervals, frequency binning). If a signal is 
detected (the threshold will in effect be more stringent than 
if n is large), the actual significance has to be estimated 
through simulations (see section 3.2). 
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